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In an alternative mechanism recently proposed, adiabatic cosmological perturbations are gener- 
ated at the decay of the inflaton field due to small fluctuations of its coupling to matter. This happens 
whenever the coupling is governed by the vacuum expectation value of another field, which acquires 
fluctuations during inflation. We discuss generalization and various possible implementations of this 
mechansim, and present some specific particle physics examples. In many cases the second field can 
start oscillating before perturbations are imprinted, or survive long enough so to dominate over the 
decay products of the inflaton. The primordial perturbations will then be modified accordingly in 
each case. 
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O ' I. INTRODUCTION 

<n : 

The detection of the cosmic microwave background (CMB) anisotropies indicates the presence of coherence over 
^ZZ^ ' super horizon scales, which is a strong indication of an early inflationary stage [1]. To date measurements are 
in agreement with the simplest prediction from inflation which is a nearly scale invariant spectrum of gaussian 
and adiabatic primordial perturbations [2]. Nevertheless, possible experimental or theoretical deviations from this 
most economical possibility have been the subject of intense study. Indeed, the increasing precision of the CMB 
measurements gives us the hope to raise the bar in the near future from the present confirmation of the general idea 
^ ' of inflation, to the study and determination of the underlying theory. 

An important aspect of inflation is the stage of reheating [3], which describes all the processes from the end of 
the inflationary expansion to the following hot "big-bang" evolution. The inflaton decay is only the first, although 
most relevant, stage of this process. The inflaton decay is typically very quick and non-perturbative, and it occurs 
immediately after the end of inflation [4]. Due to its efficiency, it leads to a distribution of particles which is very 
far from thermal equilibrium, so that reheating completes over a much longer timescale. Particle physics plays a 
. relevant role in this period. It seems conceivable that the baryon asymmetry of the universe is generated at this stage. 
On the other hand, nucleosynthesis poses strong bounds on the production of gravitationally decaying relics. These 
requirements can be combined to constrain different models of reheating. For example, limits from nucleosynthesis 
are in contrast with thermal grand unified theory (GUT) baryogenesis and only marginally compatible with thermal 
leptogenesis, while a non-thermal origin [5] can be more easily accounted for. However, while these considerations 
i-^ ■ allow us to study the inflaton interactions for any given particle physics model, we are not yet in a stage where we can 
discriminate among different scenarios. As in many other areas of physics, we are more in the need of experimental 
evidence rather than theoretical models. It is therefore important to ask whether reheating can have some other 
observational consequences which can guide us discriminating among the different possibilities. 

A positive answer to this question is provided by the recent observation that reheating could have played a key 
role also for generating the primordial perturbations leaving their imprints on the CMB, as well as on the matter 
power spectrum [6,7]. This statement challenges the common assumption that the microphysics responsible for the 
decay of the inflaton, and the successive thcrmalization, cannot have an impact on the much larger scales relevant for 
the present observations which at that time are well outside the horizon. However, perturbations generated during 
inflation can be strongly affected by the following history of the universe. For example, the adiabatic mode of the 
perturbations is sensitive to the evolution of dominating species, and it is precisely during reheating that this evolution 
is mostly unknown. Although we can be sure of the starting (unstable inflaton) and final (thermal bath) points, we do 
not know what dominates during the intermediate stages and how its equation of state evolves. Several possibilities 
can occur, depending on the inflaton potential near its minimum (for example, m 2 </> 2 versus A</> 4 potential), on the 
possible generation of unstable heavy particles which can dominate for some time, on the possible phase transitions 
and on the presence of large effective masses before thermalization completes. 

These processes can strongly affect primordial perturbations if they occur at (slightly) different moments in the 
different parts of the universe we presently observe. This can happen, for example, for the decay of the inflaton, if 
its coupling to matter is controlled by the vacuum expectation value (VEV) of some scalar field X which acquires 
inflationary fluctuations [6,7]. The equation of state of the inflaton is in general different from the one of its decay 
products, and fluctuations in the decay rate will then give rise to fluctuations in the energy density and in the 
metric. Perturbations generated by this mechanism have been named modulated perturbations. However, modulated 
perturbations do not need to be generated precisely at the inflaton decay. They could also arise due to fluctuations 
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in the decay rate of some other intermediate particle, or in the rate of the processes which are more relevant for 
thcrmalization. In general, we will denote by t* the time at which this process completes, with the assumption that 
it can be anytime during reheating. This process changes the equation of state of the universe from an unknown one 
p = wp to p = p/3 , which is typical of radiation. Due to the change in w, fluctuations in i* give rise to adiabatic 
perturbations. 

A further extension on which we want to focus, noted briefly in [7]- [9], is the possibility that the field A is evolving 
before or after t* . While A is initially frozen at some (nearly) constant value, it will start oscillating around the 
minimum of its potential (with a slowly decreasing amplitude) as soon as the Hubble expansion rate drops below its 
mass. This can significantly affect (typically decrease) the amplitude of modulated perturbations if oscillations of X 
start before perturbations are generated. Moreover, as any other modulus, X can eventually dominate the energy 
density of the universe. Isocurvature fluctuations in X will in this case be converted to curvature perturbations in 
the same way as the curvaton scenario [10,11]. Then perturbations in the new (dominating) thermal bath will replace 
modulated perturbations. A particular identification of X, and requirements for a successful reheating, may favour 
one of these cases. 

The analysis of Section II covers all these different possibilities. From the general discussion, we then special- 
ize to some specific examples motivated by particle physics in Section III. The main ingredient of the modulated 
perturbations mechansim is the existence of a scalar field with a mass less than the Hubble expansion rate during 
inflation. Theories based on supersymmerty provide the best framework for implementing this mechansim as they 
contain many scalar fields whose mass is protected against quantum corrections. We first discuss the possibility that 
A is a modulus field, with a mass ~ TeV and interactions of gravitational strength with the observable sector. This 
is probably the most natural implementation of the idea of modulated perturbations: moduli of supcrstring models 
have these properties, and their VEV's control the value of the parameters of the theory. The difficulty with this 
implementation, as we shall see, is that it typically results in a too small amplitude for the primordial perturbations, 
unless A dominates. In this case, however, one has to make sure that the late (gravitational) decay of A occurs before 
primordial nucleosynthesis. Alternative examples that we consider are the identification of A with right-handed (RH) 
sneutrinos and with flat directions of the minimal supersymmetric standard model (MSSM) , where the decay of A 
occurs at an earlier stage. In general, other scalar fields (in particular those coupled to A) should have negligible 
fluctuations in order not to affect perturbations. In Section IV we discuss possible ways to achieve such a suppres- 
sion leading us to some "benchmark scenarios" of modulated perturbations. The results will be summarized in the 
concluding Section V. 

II. GENERAL ANALYSIS 

In this section we discuss how modulated perturbations are sensitive to different assumptions on the dynamics of 
A . 1 As discussed in the Introduction, we assume that the VEV of A controls the rate r* = H (t*) of some process 
at reheating. We can in general expect 

r, = r + A"/(A) , (i) 

where T represents any A-independent contribution to that process. Here n is an 0(1) real number, while / is 
cither a (dimensionful) constant or a function which depends only weakly on A. Eq. (1) holds for a number of 
possibilities. For example, it can be applied when A is included in the vertex representing some interaction, or when 
A determines the mass of a decaying particle. For specific examples, consider the following cases of perturbative decay 
and preheating. In a perturbative decay, T* is a linear function of mass and a quadratic function of coupling. Then 
linear dependence of the coupling and mass on A will result in n = 2 and n = 1 respectively, and / is a constant. 
Whilst, for preheating to bosons, 2 T* is a logarithmic function of coupling and linear function of mass [4]. Therefore, 
linear dependence of mass and coupling on A will result in re = 1 (/ being a constant) and n = (/ a logarithmic 
function) respectively. 

Modulated perturbations are generated at t* with an amplitude 



1 A simliar analysis when X is the curvaton field has been performed in [12], [13]. 

2 The situation will be different for fermionic preheating [14], [15]. However, in supersymmetry (which provides the natural 
framework for implementing the modulated perturbations mechanism) preheating to fermions is always accompanied by bosonic 
preheating which is much more efficient. 
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The overall proportionality factor can be relevant, but we prefer to leave it unspecified not to affect the generality of 
the discussion. We assume that the potential for X has a minimum at X Q , and that only the first (quadratic) term 
in an expansion series around Xq is relevant: 
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where mx is the mass of X. We shall notice that there can be some degeneracy between Xo and To . For example, 
consider the case where perturbations are generated at the inflaton decay 3 and there are two decay modes, one 
mediated by X and one X— independent, with the same final state: C D X(jripip/M + htfiipip (ip being a fermion). 
The redefinition X = hM + X then sets To = , while changing the value of Xo . In what follows, To stands for the 
X-independent part of which cannot be removed by a redefinition of Xo. This happens, for example, when the 
two decay modes of the inflaton have different final states. 

Let us denote the Hubble rate during inflation by Hj. For mx > Hi , the field X settels at the minimum Xo with 
negligible fluctuations. We then recover the standard situation of a constant rate r = T (X ) . In the opposite case, 
mx < Hj and fluctuations of X accumulate in a coherent expectation value X/, with a dispersion SX ~ Hi . This 
gives a typical result for modulated perturbations [6,7]: 



C ~ Hj/Xi 



(4) 



As discussed in the Introduction, there are however interesting cases in which, due to the evolution of X prior or 
after t* , eq. (4) does not apply. The amplitude of modulated perturbations can be easily computed case by case and, 
although there are several parameters in the model, the final results can be most effectively summarized as a function 
of the mass and the decay rate of the scalar field X . The two cases Xo 3> X/ — Xo| and Xo <C \Xj — Xo , see fig. (1), 
lead to different results, and hence we will consider them separately in the next two subsections. 



3 Henceforth, we occassionally mention the inflaton decay as an explicit example. Nevertheless, the main conclusions will hold 
for other processes generating modulated perturbations. 
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FIG. 2. The amplitude of perturbations when Xo 3> \Xi — Xo \ . In cases (a) and (b) perturbations are generated before and 
after the start of X oscillations respectively. Case (c) corresponds to X-dominance which is essentially the same as the curvaton 
scenario. The boundary of (c) with (a) and (b) is defined by Fx = H eq , with H eq given in eqs. (5) and (6) respectively. 



A. Case I: \Xi - X \ < X < M p 

Figure (2) summarizes the results for this case. We always demand that Tx < F*, with Tx being the decay rate 
of X, which is a necessary assumption for the generation of modulated perturbations. 4 We also assume that the 
X-dependent term dominates T* and Tq can be neglected in eq. (1). For mx < F» , the field X is frozen at a value 
X/ ± Hi when the modulated perturbations are generated. The assumption Xo < M p ensures that its energy density 
is subdominant at this time. 5 This will continue to be the case provided X decays sufficiently quickly. The modulated 
perturbations are then estimated as £ ~ Hj/Xq . However, if X survives long enough, its energy density can overcome 
that of the inflaton decay products which, as discussed in the Introduction, redshifts as radiation for t > i* . While 
X is initially frozen, as soon as the Hubble parameter drops below mx , the field starts oscillating about Xo . The 
amplitude of the oscillations decreases as a -3 / 2 , where a is the scale factor of the universe, so that the energy density 
of X redshifts as that of matter. The X-dominance occurs when the Hubble parameter is 

/ Xj — X \ 4 

H cq = m x y— — -J , m x < T* (5) 

If Tx < H cq , the field X dominates before decaying. As a consequence, the relevant primordial perturbations come 
from fluctuations in the (potential) energy density of X ; px oc (Xj — Xo ± Hi) 2 , and hence £ ~ H*/\Xi — Xq| in 
this case. 

If mx > , the analysis is complicated by the fact that the oscillations of the X field start before . Dominance 
of X now occurs at 

l-3u 



— )^ rn x ( Xi ' X q ^ 4 
m x J \ M p 



Hc °l = ( ) mx ( 77 ) ' mx > r * ( 6 ) 



4 If X decays very fast and its decay products quickly thermalize, temperature of the resulting (subdominant) thermal bath 
will have superhorizon fluctuations. This will lead to fluctuations in thermal corrections to the masses and couplings of the 
fields which are coupled to the thermal bath. Modulated perturbations can still be generated if the relevant process is controlled 
by some of these fields, but we do not consider this possibility here. 

5 Mp = 2.4 x 10 18 GeV is the reduced Planck mass. 
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FIG. 3. The same as fig. (2) but for Xq <C |Xj — Xo| . Here F*j and T*^ denote the initial value of F, and its value in the 
vacuum respectively. In case (b), domination of F,,o by X-dependent and X-independent terms, see eq. (1), results in b — 1 
and b — n respectively. 



where w is the (unknown) equation of state of the inflaton decay products for r* < H < mx ■ If the equation of state 
of the universe changes throughout this interval, the first term on the right-hand side of (6) will be replaced by the 
product of terms expressing the redshift in each sub-interval. If Tx < H cq , the X field dominates and perturbations 
in the energy density of X will give again £ ~ Hi/\Xi — Xq \ . For Tx > H eq , however, modulated perturbations 
will have a smaller amplitude than the one found above due to the fact that 5X also decreases as a -3 / 2 during the 
oscillations of X . 

To summarize, the amplitude of primordial perturbations will amount to 

f (-)* < f < v^x-v < 7 » 

X \m x J X \Xi-X \ 

in cases (b), (a) and (c) respectively. 6 Since \Xj — Xq\ <C Xq, obtaining acceptable perturbations requires that 
Hj/\Xj — Xo\ > 10 -5 . On the other hand, current limits on the non-gaussianity of perturbations set the upper bound 
Hj/X Q < 1CP 1 [16]. Note that for values close to this limit, perturbations of the correct size can be generated only 
when r* < mx- 

These results should be taken only as estimates of the exact value. More precise values can be obtained once the 
rate (1) is specified. Clearly, one should not expect a discontinuity of The amplitude of the adiabatic mode will 
smoothly interpolate between the values given in fig. (2) and the bordering regions are characterized by a significant 
amount of isocurvature perturbations. 



B. Case II: X < \Xi - X \ < M p 

We now discuss the case where the initial displacement of X is much greater than X . In addition, we focus on the 
situation in which the rate (1) is initially much greater than in the vacuum of the theory, i.e. that r*,j = T* (Xi) 3> 
r^o = r* (-^o) • L e t us first discuss the simpler case mx < I\/ , characterized by a frozen field X at i» . In this 



6 If Fq dominates in eq. (1), perturbations will be further suppressed by a factor of Xq f(Xo)/Vo in cases (a) and (b), while 
remaining unchanged in case (c). 
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case, the analysis is analogous to the one in the previous subsection. If the field A never dominates, the amplitude of 
modulated perturbations now amounts to £ ~ Hj/ ' Xj . Remarkably, the same value (up to an 0(1) factor) is obtained 
also if the field X survives long enough to dominate, since p ~ px oc (Aj ± Hi) 2 in this case. 

The situation for mx > F*,i is instead more complicated. As soon as X starts oscillating with a decreasing 
amplitude, the rate T* changes according to 

T*j ^T(X)ocX n <xH^ , (8) 

where we remind that w denotes the equation of state before t* . If n > 1 + w (expected, for example, for perturbative 
processes) T» decreases more rapidly than (or at the same rate as) H . Hence, the process which generates modulated 
perturbations will not be efficient until its rate stabilizes at r».o. Note that the situation in this case is very different 
from the one discussed in the previous subsection, where r* ~ r*^ at all times. At t = i* , the field X evaluates to 

X(U)~(X I ±H I )(^) 1+a . (9) 
\m x J 

First assume that r* ;0 ~ Xff /(A ), i.e. that the A-dependent term dominates in the vacuum. Since f{X) is assumed 
to be a slowly varying fucntion of X, see the discussion after eq. (1), Xi/Xq ~ (r».//r».o) 1//l1 when n > 1 + w. Thus 
the expression for the ampltiude of perturbations can be cast in the form 

TT / T \ l/ n / | 1 

til I J- *,/ \ / 1 *.0 



- (10) 

The situation will be somewhat different when the A-indcpendcnt term dominates in the vacuum. Now r* ; o ~ To, 
while <5r* o ~ A} l " 1 i7 / (r* ! o/™x)" /1+tu . This results in a smaller value for the amplitude of perturbations 

Hi /r*,/ \ i r*,o 



If A dominates, we have again £ ~ Hi/Xi . This occurs provided A decays after H eq , given by eq. (6). 

It turns out from (10) and (11) that ( < Hi/Xi in case (b). It is therefore possible to generate acceptable 
perturbations in this case for 1CP 5 < Hi/Xi < 10 _1 . The upper bound is again given by current limits on the 
non-gaussianity of perturbations. When Hi/Xi is close to this value, the modulated perturbations mechanism can 
work even if A/ is entirely due to the accumulation of quantum fluctuations, and if inflation does not last much longer 
than 60 e-foldings [17]. Note that the last two terms on the right-hand side of (10) and (11) should yield a number 
> 10~ 4 . Otherwise oscillations of A will suppress modulated perturbations to an unacceptably low value. 

III. PARTICLE PHYSICS EXAMPLES OF X 

In this section we present three realistic particle physics examples of A. We will discuss different possibilties for 
generating sufficient perturbations while satisfying other cosmological constraints, particularly to avoid the gravitino 
problem. As we will point out, these examples actually represent a wide range of particle physics candidates for A. 

A. Moduli 

In superstring models physical parameters are usually set by the VEV of moduli scalar fields, which acquire mass 
only after supersymmetry is broken. In a cosmological setting, supergravity effects typically provide these fields with 
a mass of the order of the Hubble parameter H [18,19]. However, there are situations in which (for example, due to 
a "Heisenberg symmetry" [20]) supergravity corrections do not affect some of these directions. These fields may then 
naturally be expected to have a mass of the order of the low energy supersymmetry breaking; mx ~ O(TeV). Moduli 
are gravitationally coupled to the other fields in the theory, so that their decay rate is estimated to be Tx ~ m x /Mp. 
They are typically characterized by A ~ M p , thus leading to case I in above. 

Modulated perturbations can be expected if any parameter relevant for reheating is controlled by one of these 
fields [7]. However, as we shall see, they turn out to be rather small in this example. The final result for the 
perturbations is sensitive to the difference |A/ — A | , namely on the initial displacement of the modulus from the 
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minimum of the potential. We regard this quantity as a free parameter, with the only constraint that it should be 
greater than the Hubble parameter during inflation Hi. From eq. (5) we see that for \Xj — Xq\ < 10 11 GeV the 
modulus field decays before it dominates. However, as we just noted, the limit on \Xo — Xj\ povides an upper bound 
also on Hj. This results in modulated perturbations that are too small, see fig. (2). 

For \Xj — V |>10 1:L the modulus decays after it dominates, and the relevant perturbations are the ones in the energy 
density of X . 7 As long as \Xi — Xq\ < M p , the amplitude of the perturbations amounts to ~ Hi/\Xj — Xq\ (see 
fig. 2). Acceptable perturbations are then generated provided Hi ~ 10~ 5 \Xi — Xq\ . Actually mx ^ 50 TeV is needed 
so that the moduli decay will result in a reheat temperature ~ O(MeV), compatible with the big bang nucleosynthesis. 
This can be achieved while keeping the mass splitting between matter fields and their supersymmetric partners at 
TeV level, for example, in models of anomaly-mediated supersymmetry breaking [21]. For |Aj — X \ > M p , the 
modulus X drives a stage of inflation, and the analysis of the previous section does not apply. However, due to the 
smallncss of mx,i.e. mx <C 10 13 GeV, identification of modulus with the inflaton field results in unacceptably small 
perturbations. 



B. Right-handed sneutrinos 

The right-handed (RH) sneutrinos arise in supersymmetric extensions of the standard model [22] which explain the 
smallness of the mass of the left-handed (LH) neutrinos via the see-saw mechanism [23]. A non-zero VEV for N at 
the minimum of its potential breaks -R-parity, thus destabilizing the lightest supersymmetric particle (LSP) [22]. This 
VEV should be ^ 1 GeV, if the LSP survives until today and constitutes the dark matter. Nevertheless, since TV is a 
standard model gauge singlet, it is conceivable that the V-indcpdcndent contribution to the process which generates 
modulated perturbations can be recasted in a much larger X , as we have discussed after eq. (3). Consider an example 
where perturbations arc imprinted at the inflaton decay which proceeds through the superpotential couplings (the 
boldface charachters denote superfields) 

+ (12) 

M p 

Then Xq = hM pi leading to case I (II) in above if Xi <C hM p (^ hM p ). We concentrate on the latter possibility as 
the former will lead to a situation similar to that discussed in the previous subsection. 8 
The see-saw formula gives [23] 

m x > y 2 (H) 2 /m„~y 2 l0 15 GeV , (13) 

where (H) ~ 200 GeV is the electroweak symmetry breaking scale, m v ~ 0.05 eV is the mass of the heaviest light 
neutrino and y denotes the neutrino Yukawa coupling. The same interactions give rise to a sneutrino decay rate 

Fx = z— m x ~ J/ 4 10 14 GeV . (14) 

In realistic models of neutrino masses based on the see-saw mechansim [24], mx is typically much larger than the 
scale of electroweak symmetry breaking. Moreover, producing sufficient baryon asymmetry via leptogenesis [25] leads 
to additional constraints on the model parametrs. In particular, various leptogenesis scenarios [26-30] require that 
nix > 10 5 GeV (for thermal leptogenesis the bound is about 4 order of magnitude stronger) unless some specific 
fine-tuning occurs (for example, having nearly degenerate RH (s)neutrinos [31]). Therefore, nix 3> 1 TeV in general. 

In principle, T* can be larger or smaller than mx ■ In the former case, modulated perturbations must be generated 
by a very rapid process, most notably non-perturbative inflaton decay via preheating. A large may be associated 
with a high reheat temperature Tr , 9 which can then lead to thermal overproduction of gravitinos. In gravity- 
mediated models of supersymmetry breaking the gravitino mass is m 3 / 2 — 100 GeV — 1 TeV and (up to logarithmic 
factors) [33] 



7 To be precise, nix ^ 50 TeV is needed so that the moduli decay will result in a reheat temperature > O(MeV), compatible 
with the big bang nucleosynthesis. This can be achieved while keeping the mass splitting between matter fields and their 
supersymmetric partners at the TeV level, for example, in models of anomaly-mediated supersymmetry breaking [21]. 

8 Generating density perturbations from the RH sneutrinos has also been considered in [32]. 

9 We have defined i„ as the moment when the inflaton decay products have an equation of state of radiation. This does not 
necessarily mean that their distribution is thermal. Once thermalization of decay products has completed, Tr will be the 
highest temperature in the radiaiton-dominated era. 
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Bounds from nucleosynthesis (due to photodissociation of light elements by gravitino decay products) impose T R < 10 9 
GeV [34], corresponding to Hr < 10 GeV (Hr being the Hubble parameter at T = Tr). In models of gauge-mediated 
supersymmetry breaking, the gravitino is the LSP and m 3 / 2 can be as low as 1 KcV. Its fractional energy density is 
in this case [35,36] 

where h w 0.7 and M3 is the gluino mass parameter. Then, for a typical value M3 ~ 500 GeV, the dark matter 
limit results in Tr < 10 s m 3 / 2 , corresponding to iJ^ < 10 17 m^ 2 /M p . In both cases higher Tr can be accommodated, 
provided the gravitino abundance is diluted by an entropy release at later times. This entropy can be naturally 
provided by the decay of the sneutrino if its energy density becomes significantly greater than that of the thermal 
bath. In addition, the observed baryon asymmetry of the universe can also be generated by such a late decay [29,30]. 
Then (5), (14), (15) and (16) lead to the bounds 



»<5xio»r^v^2| 



and 



M p 

in gravity-mediated and gauge-mediated models respectively. We are thus naturally led to consider the case in which 
the sneutrino dominates. Primordial perturbations will then be related to fluctuations in the energy density of X , 
corresponding to case (c) in fig. (3). 

Let us finally consider the case where r*,/ < mx ■ If the reheat temperature respects the gravitino bound, no entropy 
generation will be required. Then the sneutrino may decay before or after dominating the universe corresponding to 
cases (b) and (c) in fig. (3), respectively, and perturbations will be given by (10) in the former case. If the reheat 
temperature exceeds the gravitino bound, sneutrino dominance will be favoured again to solve the gravitino problem. 



C. Supersymmetric flat directions 



There are many directions in the space of the Higgs, slepton and squark fields along which the scalar potential 
of the MSSM identically vanishes in the limit of exact supersymmetry [37]. These flat directions acquire a mass 
from supersymmetry breaking. In gravity- mediated models m x = 100 GeV — 1 TcV [22]. In gauge- mediated models 
mx — 100 GeV — 1 TeV for small (X), it drops oc (X)^ 1 at intermediate VEVs and mx ~ m 3 / 2 for large (X) [36]. 

The Higgs fields have a VEV ~ 100 GeV in the vacuum, while those of the sleptons and squarks vansih. Therefore 
identification of X with the MSSM flat directions (the primary example considered in [6]) leads to case II in above 
where perturbations read from fig. (3). The situation is qualitatively similar to the sneutrino example but some 
differences exist. For the flat directions mx is typically much smaller than the RH sneutrinos, particularly in gauge- 
mediated models. This allows a smaller T* (as expected for perturbative processes) when r*,j > mx- The fact that 
the Higgs, slepton and squark fields are gauge non-singlets can also result in a different situation when IV / < mx- 
Consider again the example of inflaton decay in eq. (12), with X being a Higgs, slepton or squark field. Now the 
two terms in (12) must couple the inflaton to different final states and ^2, resulting in decay rates Tx and Tq 
respectively. In any acceptable scenario r„ i0 — T , and hence perturbations will be given by (11), instead of (10) for 
the sneutrinos, provided X does not dominate. 

The V-dominance leads to case (c) in fig. (3). The decay of the flat directions proceeds through the SM Yukawa 
and gauge couplings. In particular, one may wonder whether X can dominate at all as gauge couplings lead to a 
rapid decay of X. However, the gauge (and Yukawa) couplings also result in an effective mass m e ff for the decay 
prodcuts due to their couplings to the X condensate and thermal effects. The decay is kinematically forbidden so long 
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as m e ff > mx [39]. It occurs when H = = mhi[Tx, Hki n ], with Hki n denoting the value of Hubble parameter at 
which m e ff drops below mx- 10 Therefore A-dominance actually requires that Hd < H eq , where H eq is given by (6). 

D. Summary 

In supcrgravity models all scalar fields typically receive a supersymmetry breaking mass ~ 7713/2. If a supcrsym- 
metry conserving (supcrpotential) mass term is allowed for a field, its mass can be 3> 7713/2- However, a tree-level 
(supersymmetry breaking or conserving) mass term will be forbidden if the theory has some symmetry in the super- 
potential (thus resulting in a Goldstone boson) or in the kinetic function (as in no-scale supergravity [40]). In this 
case the corresponding field obtains a mass <C m 3 / 2 through symmetry breaking interactions at higher orders. Thus, 
unless forbidden by some symmetry, mx > 777,3/2 will be expected in general. Under the most general circumstances, 
a massive scalar field is either in a hidden sector, thus gravitationally coupled to the SM fields, or has gauge and/or 
Yukawa couplings to the observable sector. Moduli are an example of the case with mx ~ m 3 / 2 an d gravitational 
decay to matter fields. The supersymmetric flat directions (mx > 777,3/2) and the RH sneutrinos (mx 3> m + 3/2) 
have gauge and Yukawa couplings to the SM fields. The three examples considered above therefore represent a wide 
range of possible X candidates. 

One comment is in order before closing this section. We have only considered the mass term in V(X), see (3). 
However, higher order terms are naturally expected to arise and make the potential steeper than X 2 at large field 
values. This will result in an upper limit X max on Xj at which V"(X max ) ~ Hj. When Xj ~ X maxi fluctuations of 
X are attenuated even before oscillations start [41]. This will suppress the amplitude of perturbations compared to 
that is given in figs. (2) and (3), which will be valid so long as Xi < X max . A limit on Xi, through the expressions 
for the amplitude of perturbations, translates into an upper bound on the scale of inflation Hi. The strongest bound 
Hi < 10~ 5 X max is obtained when ( ~ H1/X1, while cases with £ < H1/X1 lead to weaker constraints on Hi. In 
particular, Hi can be just one order of magnitude below X max when the value of H1/X1 is close to the current 
limit from non-gaussianity of perturbations. The bound on Hi can in this case be comparable with (or weaker 
than) Hi < 10 13 GeV, which is required so that the inflaton fluctuations not yield too large perturbations, even if 
X m ax *C M p . The value of X maxi signifying the importance of higher order terms, strongly depends on the particle 
physics identification of X. For a modulus field these terms are M p suppressed and the mass term will be dominant 
up to Xi < M p . For the RH sneutrinos and supersymmetric flat directions, X max can assume a wide range of values 
<C M p . Due to the model-dependence of X max and the constraint it imposes on Hi, we have not considered higher 
order terms in V(X) in our discussion. 

IV. DISCUSSION 

So far we have assumed that only the inflaton and the X field dynamically evolve in the early universe. Any other 
scalar fields with masses smaller than Hi are also expected to acquire inflationary fluctuations, and be substantially 
displaced from the minimum of their potential at early times. They can then contribute to density perturbations in 
a similar fashion as the X field. Since X is considered to be the dominant source for generating perturbations here, 
fluctuations of other scalar fields should be suppressed. This is particularly true for X decay products, generally 
denoted by x- This can be naturally achieved if these fields have a mass >Hi during inflation. One possibility is that 
supergravity effects provide such a mass [18,19]. Note, however, that the supergravity mass for X should be < Hi. 
As pointed out earlier, this can happen if X has a different kinetic function from other fields. It is also possible that 
different masses for different fields arise dynamically as a result of quantum corrections [42,43]. Fluctuations of x will 
also be suppressed if its coupling to X and/or inflaton, denoted by y and h respectively, is sufficiently large, i.e. that 
yXi > Hi and h(f>o > Hi respectively (</>o is the initial amplitude of inflaton oscillations). Since mx < Hi, preheating 
decay of X and/or inflaton to x will be expected in this case [4]. 

The discussion of Section III, together with the requirement for suppressing fluctuations of x, lead us to the following 
"benchmark scenarios" for modulated perturbations: 



For many flat directions oscillations of the zero-mode condensate can be fragmented into non-topological solitons, so called 
supersymmetric Q-balls [38], which may be long-lived or stable. This will further complicate the situation, and hence we do 
not consider it here. 
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1. X is in the hidden sector and mx ^ 50 TeV. Supergravity effects provide a mass > H for \ suppressing its 
fluctuations. X dominates the universe and fluctuations in its energy density give rise to perturbations; case (c) 
in fig. (2). The situation is essentially the same as the curvaton scenario. 

2. X is in the observable sector and mx ^ 1 TeV. Perturbations are generated when the inflaton decyas via 
preheating, before X starts oscillating; case (a) in fig. (3). Coupling of X to \ is large enough (y > 10~ 5 ) to 
suppress its fluctuations. This also leads to a rapid non-perturbative decay of X to \i which is welcome as X 
will not dominate and consequently affect perturbations. The universe thermalizes sufficiently late, thus there 
will be no thermal overproduction of gravitinos and no late time entropy generation will be required. 

3. X is in the observable sector and starts oscillating before perturbations are imprinted; case (b) in fig. (3). This 
happens due to the slowness of the process generating perturbations, for example, gravitational decay of the 
inflaton or decay of a long-lived massive particle produced in the inflaton decay. Such a slow process will be 
required by the gravitino bound when the universe thermalizes rapidly. Oscillations of X in this case suppress 
modulated perturbations compared to that in the previous scenario. Preheating decay of X is not allowed, thus 
yXi < mx, since it should survive long enough until perturbations are generated. This rules out suppressing 
X fluctuations through its coupling to X (supergravity effects or a large coupling to the inflaton can do this 
instead). On the other hand, since fluctuations in the energy density of X can be 3> 10~ 5 , X must decay long 
before dominating in order not to affect perturbations. This, after using (6), will result in an upper bound on y. 

We note that X can undergo early oscillating if it acquires a large thermal [44] or non-thermal [45] mass mx, e ff- 
This will happen provided \ are m equilibrium with the (p)reheat plasma and y is sufficiently large, 11 so that 
m x,eff exceeds the Hubble parameter when H > mx- 

4. X is in the observable sector and starts oscillating before or after modulated perurbations are generated. The 
reheat temperature exceeds the gravitino bound but X dominates the universe before decaying, thus diluting 
the gravitinos in excess. Fluctuations of \ wm m this case be suppressed by supergravity effects or a large 
coupling to the inflaton. Perturbations will be due to fluctuations in the energy density of X; case (c) in fig. (3). 

V. SUMMARY AND CONCLUSIONS 

In this paper we have studied different scenarios of modulated perturbations. In this mechanism adiabatic cosmo- 
logical perturbations are generated during reheating when the equation of state in the different parts of the universe 
changes at different moments. The inhomogeneity arises due to inflationary fluctuations of some scalar field X whose 
VEV controls the value of parameter(s) (such as mass or coupling) involved in the relevant process. This is a generic 
property expected in models based on superstring and supersymmetric theories. Depending on the values of the rate 
of the process T*, mass of the scalar field mx and its decay rate Fx, different situations can arise. We presented a 
general analysis of possible scenarios, with the main results summarized in fig. (2) and fig. (3). We also introduced 
some specific particle physics examples including the string moduli, RH sneutrinos and supersymmetric flat directions. 
These examples, in which mx is of the same order or greater than the gravitino mass 1713/2, represent a wide range 
of particle physics candidates of X. 

The simplest possibility is that perturbations are generated when X is frozen at an initial value X\. If mx > 
O(TeV) (as for the RH sneutrinos and supersymmetric flat directions), this requires a very rapid process, namely 
non-perturbative inflaton decay via preheating. Moreover, X must decay before dominating the energy density of the 
universe, and the reheat temperature should be sufficiently low to avoid overproduction of gravitinos. For a slower 
(likely perturbative) process, X can start oscillating around the minimum of its potential before perturbations are 
imprinted. Then the expansion of the universe redshifts fluctuations of X and dampens perturbations. An interesting 
consequence is that the level of non-gaussianity can in this case be just below the current limits and accessible to 
future experiments. If X is long-lived (for example, a modulus field), it can dominate the universe before decaying. 
In this case the situation will essentially be the same as the curvaton scenario and fluctuations in the energy density 
of X give rise to perturbations. The X-dominance can be a virtue when the reheat temperature exceeds the gravitino 
bound necessitating late time entropy production. The amplitude of the adiabatic mode is in general very different 
from one case to another. It interpolates between the different values and bordering regions are charachterized by a 
significant isocurvature component. 



y should not be too large, however, otherwise re-scatterings by \ quanta will quickly destroy the X condensate. 
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In conclusion, the modulated perturbations mechanism is a viable alternative with potentially interesting observa- 
tional consequences. It can be implemented in various ways and a successful particle physics identification can lead 
us to a clearer and more complete picture of reheating. 
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